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Abstract 

A link between memory effects in quantum kinetic equations and nonequi- 
librium correlations associated with the energy conservation is investigated. 
In order that the energy be conserved by an approximate collision integral, 
the one-particle distribution function and the mean interaction energy are 
treated as independent nonequilibrium state parameters. The density opera- 
tor method is used to derive a kinetic equation in second-order non-Markovian 
Born approximation and an evolution equation for the nonequilibrium quasi- 
temperature which is thermodynamically conjugated to the mean interaction 
energy. The kinetic equation contains a correlation contribution which ex- 
actly cancels the collision term in thermal equilibrium and ensures the energy 
conservation in nonequilibrium states. Explicit expressions for the entropy 
production in the non-Markovian regime and the time-dependent correlation 
energy are obtained. 
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I. INTRODUCTION 



Experimental studies of fast relaxation processes caused by the interaction of short laser 
pulses with matter have inspired a renewed interest in non-Markovian kinetics. Al- 
though this subject has been under development for many years, recent investigations have 
shown that the inclusion of memory effects in collision integrals leads to some serious prob- 
lems, which did not receive proper attention previously. 

Since in many cases of experimental interest the system can be described in terms of 
weakly interacting quasiparticles, the most-used non-Markovian quantum kinetic equations 
are in fact modifications of the so-called Levinson equation [H in which the collision integral 
is taken in second-order Born approximation and the energy-conserving delta function is re- 
placed by an oscillating memory kernel (see, e.g. 0). It should be noted, however, that such 
kinetic equations have some grave disadvantages, (i) The Levinson-like kinetic equations 
have unstable solutions and even produce negative distribution functions, (ii) If the initial 
state of the system is already thermal equilibrium, the collision integrals do not vanish, 
giving rise to non-physical time evolution. The first defect can be overcome by using certain 
decaying memory kernels which take account of the quasiparticle damping. However, 

in such approaches, the problem of the equilibrium solution still persists because all "im- 
proved" memory kernels lead to collision integrals which do not conserve the total energy 
of the system. As a consequence, in the long-time (Boltzmann) limit the "improved" non- 
Markovian quantum kinetic equations do not describe the relaxation to thermal equilibrium, 
leading to an overpopulation of high-energy states. 

To summarize, it appears that the quasiparticle damping alone cannot be responsible for 
the long-time asymptotic behavior of non-Markovian quantum kinetic equations; there must 
be another physical mechanism which provides the exact cancellation of collision effects in 
thermal equilibrium and does not violate the energy conservation. The origin of such a 
mechanism can be deduced from consideration of kinetic processes in the presence of initial 
correlations. As early as 1970 Lee et al. studied the evolution of a weakly interacting 
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low-density classical gas with a correlated initial state and showed that in equilibrium the 
changes in the one-particle distribution function due to collisions (including memory effects) 
and initial correlations exactly cancel each other. Recently the same result was obtained for 
quantum systems 0]. It therefore would appear reasonable that the interplay between col- 
lisions and nonequilihrium correlations would be the mechanism which provides the correct 
long-time behavior of non-Markovian quantum kinetic equations. 

This paper presents an approach in which the time evolution of the one-particle dis- 
tribution function is coupled with the evolution of long-lived correlations associated with 
conservation laws. The basic idea is to treat conserved quantities, most notably the energy, 
as independent state parameters in addition to the one-particle distribution function. Al- 
though this is not to say that all many-particle correlations can be incorporated in this way, 
the advantage of the above idea, originally suggested in the context of a generalization of 
the Enskog theory to dense classical gases with "soft" inter-particle potentials [H,^ and then 
applied to quantum systems pO| , p!l|] , lies in the fact that now, due to microscopic equations 
of motion, the energy is an exactly conserved quantity in any approximation for the collision 
integral. Moreover, the energy conservation leads to the appearance of additional correlation 
terms in the kinetic equation, which substantially compensate the collision contribution. In 
this paper, within second-order non-Markovian Born approximation, we derive a collision 
integral involving the collision and correlation terms. It is shown that the kinetic equation 
is consistent with the energy conservation and has the correct equilibrium solution. 

The paper is organized as follows. In Section II we consider a generalized Gibbs en- 
semble in which the energy and the one-particle distribution function play the role of given 
nonequilihrium state parameters. The corresponding relevant statistical operator is used to 
define the thermodynamic quantities: nonequilihrium entropy, the quasi-temperature, and 
the quasi-chemical potential. We also derive the evolution equations for the thermodynamic 
quantities in terms of the collision integral. Section III is concerned with construction of the 
nonequilihrium statistical operator describing non-Markovian kinetic processes in a spatially 
homogeneous system. The nonequilihrium statistical operator is found in terms of the rele- 
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vant statistical operator for quantum systems with a weak interaction by employing iterative 
solution of the von Neumann equation. In Section IV the nonequilibrium statistical operator 
is used to calculate the non-Markovian collision integral. In Section V we derive the entropy 
production in the non-Markovian regime and an expression for the time-dependent correla- 
tion energy of the system. Section VI sketches a generalization of the theory to spatially 
non-homogeneous systems. Finally, Section VII contains the final conclusions and comments 
on further applications of the theory. 

II. NONEQUILIBRIUM CORRELATIONS ASSOCIATED WITH 

CONSERVATION LAWS 

To put our discussion into a straightforward language, we shall consider a system of 
fermions or bosons with the Hamiltonian H = Hq + H', where Hq is the kinetic energy 
operator and the term H' describes a pairwise interaction between particles. In second 
quantized form, these operators are given by 

Ho = J^i^u'ahi', H' = - {I' 2' \V \12) a\,a\,a^a2, (2.1) 

1 ^ 121'2' 

where the label k denotes a complete set of single-particle quantum numbers, her- 
mitian single-particle energy matrix, a^. and a\ are Fermi or Bose annihilation and creation 
operators. Generally speaking, the Hamiltonian can also contain additional terms describ- 
ing interaction of the system with external fields. For simplicity, we will not consider this 
interaction explicitly and restrict ourselves to relaxation processes in the system just after 
the initial excitation by the external field. However, the theory can easily be generalized to 
the case where the field effects are taken into consideration (see Section VII for a discussion) . 

In what follows the system is assumed to be spatially homogeneous, which is adequate 
for most experimental situations where memory effects are of crucial importance 0. Some 
aspects of the following analysis may also be of interest for transport processes in non- 
homogeneous systems, so that we shall return to this point in Section VI. In the spatially 
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homogeneous case, it is convenient to take (k) = (p^^., Cf^), where p is the momentum and a 
is the spin index. Then the kinetic energy operator in Eq. becomes 



^o = E^i«K (2-2) 
1 

where ei = are single-particle energies. 

Within the kinetic description of the system, the main objective is to derive a kinetic 
equation for the one-particle distribution function 

h{t) = {a[a,Y ^ (2.3) 

where the symbol (■ ■ ■)* stands for the average calculated with some nonequilibrium statisti- 
cal operator Q{t). Formally, the kinetic equation follows immediately from the von Neumann 
equation for the statistical operator 

f^.i-[.W,.]^0, (.4, 

Taking the product of this equation with f-^ and then calculating the trace, we get 

dm 

dt 

where 



hit). (2.5) 



hit) = ^{[Lh']Y = ^Tt {[f„ H'] g{t)} (2.6) 

is the collision integral. In order that Eq. (|2.5| ) be a closed kinetic equation, the collision 
integral must be expressed in terms of the one-particle distribution function — which is 
equivalent to requiring that the statistical operator g{t) is represented as a functional of the 
one-particle distribution function. A common technique for obtaining such a representa- 
tion is based on the condition of complete weakening of initial correlations for the quantum 
BBGKY hierarchy in combination with some truncation procedure. A limitation of this 
standard scheme is the very strong assumption that the one-particle distribution function 
is the only nonequilibrium state variable describing the system. However, there exist long- 
lived correlations which cannot be expressed exactly in terms of the one-particle distribution 



function. For example, many-particle correlations arise due to formation of bound states. 
Another origin of long-lived correlations lies in local conservation laws. The conservation 
of energy is of special importance because the density of the interaction energy is deter- 
mined by the two-particle distribution function. Thus, strictly speaking, kinetic processes 
must always be considered together with the evolution of locally conserved quantities, i.e., 
with hydrodynamic processes. A consistent description of kinetics and hydrodynamics can 
be developed by treating the one-particle distribution function (for quantum systems, the 
Wigner function) and the energy density as independent state parameters, which means 
that the statistical operator g{t) is represented as a functional of the corresponding dynam- 
ical variables (see, e.g., |[TT|-[T3| and references therein). Here we will follow this approach 
to study non-Markovian kinetic processes in a spatially homogeneous system. In this case 
the correlation effects arise due to the fact that the total energy is an integral of motion. 
Modifications of the theory needed for the spatially non-homogeneous case will be discussed 
in Section VI. 



A. The relevant statistical operator 

We begin by considering the statistical thermodynamics description of nonequilibrium 
states with correlations caused by conservation laws. For a one-component spatially homo- 
geneous system, the conserved quantities of interest are the total number of particles, A^, 
the total momentum, P, and the total energy, S. They can be expressed as the mean values 
of the dynamical variables N = J2i fi, P = Si Pi/i) and H, i.e., 

N = {Ny = j:fiit). p = (p)* = Epi/iW> ^ = {Hy. (2.7) 
1 1 

For simplicity, we shall assume the total momentum to be zero. Then the only integrals of 
motion that must be considered are the energy and the number of particles. 

To see how the statistical thermodynamics description of nonequilibrium correlations can 
be developed, we recall the well-known grand canonical statistical operator 
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g^^ = exp {-/3{H - /iiV)} /Ti exp {-/3{H - /iiV)} , (2.8) 

where the equihbrium inverse temperature, (3 = 1/T, and the chemical potential, fi, are 
related to N and S by the equations of state 

iV = Tr {iVf., J , £ = TT{Hg^^}. (2.9) 



The grand canonical operator (|2.8| ) describes the equilibrium Gibbs ensemble and is a sta- 
tionary solution of the von Neumann equation (|2.4| ). Suppose now that we want to con- 
struct a statistical operator which describes a generalized Gibbs ensemble characterized by a 
nonequilibrium one-particle distribution function f-^(t) and by given values of the conserved 
quantities. This is a special case of a more general situation where the state of the system 
is described by a set of parameters Pm{t) which can be represented as the mean values, 
Pm{t) = Tr |-P„f)(t)|, of some dynamical variables P^. As argued by Jaynes fl^, the corre- 



sponding relevant statistical operator can be derived by maximizing the entropy functional 



for given P^if). The extremum condition for the entropy gives |]T2 



Q,,,{t)=e^v\-Y.KXt)Pm] Ar exp|-^A^(t)P„|. (2.10) 

The Lagrange multipliers A^(t) are to be expressed in terms of the P„(t) from the self- 
consistency conditions 

Pjt)=Tr{P„f.,,i(t)}, (2.11) 

which can also be interpreted as nonequilibrium equations of state. Following the above line 
of reasoning, it is easy to derive the relevant statistical operator which describes the gener- 
alized Gibbs ensemble with given values of the total energy, the total number of particles, 
and the one-particle distribution function. We write it in the form 



&el(^) = ^ 

^rel 



^— exp |-/3*(t) (h - /i*(t) n)-Y. Ai(t) A I , (2.12) 



where the partition function is determined by the normalization condition. 



7 



= Tr exp {h - ^^*{t) iv) - ^ X,{t) a| , (2.13) 

and the Lagrange multipliers P*{t), fi*{t), and A^(t) are to be calculated from the self- 
consistency conditions (nonequilibrium equations of state) 

A(t) = Tr {f,g,At)} , iV = Tr {Ng^^.it)} , £ = Tr {Hg,,,it)} . (2.14) 



Note that the grand canonical operator ( |2.8| ) is a special case of the relevant statistical 
operator ( p.l2|) , where the state parameters are integrals of motion. 



B. Alternative representations for the relevant statistical operator 

By analogy with Eq. (p.8|), the quantities T*(t) = 1/P*(t) and /i*(t) may be called the 
quasi-temperature and the quasi- chemical potential, respectively. Here one comment is in 
order. The important difference between the state parameters and S is that the former 
is a linear combination of the one-particle distribution functions, /i(t), whereas S involves 
the mean interaction energy 

S,^,{t) = {H')\ (2.15) 

which, in general, cannot be expressed in terms of fi{t) only. Thus, in fact, the set of 
independent state parameters consists of fi{t) and £-^^^^{t). This also is seen from Eq. (p. 12 



where the terms coming from the kinetic energy operator and the particle number operator 
can be combined with the last term by introducing new Lagrange multipliers h.^{t) through 
the relation 

Ai(t)=/5*(t)[ei-/i*(t)] + A,(t). (2.16) 
Then Eq. ( |2.12| ) reduces to 



QAt) = exp|-/3*(t)^'-^Ai(t)/;| (2.17) 

with the partition function given by 
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= Tr exp ^' - ^ A,{t) a| . (2.18) 

The self-consistency conditions for the Lagrange multiphers (3*{t) and A^lt) can now be 
taken in the form 

= Tr{/;ft,i(t)}, £,,,it) = TT{H'g,,,it)}, (2.19) 

since the mean kinetic energy is exactly expressed in terms of the one-particle distribution 
function. Summarizing, the quasi-chemical potential, /U*(t), can, in principle, be excluded 
from the set of Lagrange parameters since the second of the conditions ( |2.14| ) follows from 
the first one. Nevertheless, the expression (|2.12|) has two advantages. First, it goes over 



explicitly to the equilibrium statistical operator, if P*(t) = 1/T, n*{t) = /x, and \i(t) = 0. 
Second, some formulas to be derived in the following have a more clear physical interpretation 
when written in terms of the quasi-chemical potential. For these reasons, we will consider 
Eq. ( p.l2|) as a representation for Qj-^iit) which is equivalent to Eq. (|2.17| ) by virtue of 
Eq. ( CT ). Note, however, that Eq. (|2l6|) determines only X-^{t) - p* (t) /j,* (t) , but not 
the quasi-chemical potential itself. Since no physical quantity depends on the special choice 
of the quasi-chemical potential, the function fi*{t) may be chosen arbitrarily provided that 
fi*{t) = fi in thermal equilibrium. For our purposes, it will be convenient to define the 
quasi-chemical potential through the condition 

N = TT{Ng^{t)}, (2.20) 

where ^q(t) is the auxiliary quasi- equilibrium statistical operator 

0^{t) = exp {-P*it) [h - /i*(t) n)] /Tr exp {-(3*{t) [h - fx*{t) n)} , (2.21) 

which describes the state characterized by the quasi-temperature T*(t) = l/j3*{t) and the 
total number of particles A^. Equation ( |2.2CI| ) ensures that, in the equilibrium limit, fi*{t) 
goes to the chemical potential fi since l3*{t) goes to the equilibrium inverse temperature. 

For completeness, we give one more representation for the relevant statistical operator, 
which is obtained from Eq. (p.l2|) by the formal decomposition of the Hamiltonian 
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H = n^{t)+n'{t), (2.22) 

where 

n,it) = ^E,{t)a\a„ n'{t) = H'-^Ef''it)a\a„ (2.23) 
1 1 

and the re-normahzed single-particle energies, 

E,{t)=e, + Er{t), (2.24) 
involve the exchange Hartree-Fock term 

Sr^W=E(12|V^|12)e./2(t). (2.25) 

2 

Here and in what follows the subscript "ex" indicates the symmetrized (antisymmetrized) 
interaction amplitude 

(12|r|l'2')^^ = (12|V|1'2') T (12|V|2'1') = (12|V|1'2') t (21|V|1'2') (2.26) 



with the upper sign for fermions and the lower sign for bosons. Insertion of Eq. ( |2.22|) into 
Eqs. ([2A^ ) and (g7[3|) gives 



ftei(i) = exp |-/5*(t)7i'(t) - E Ai(t) /i| , (2.27) 

Z,,i(t) = Tr exp i^-(3*{t)n'{t) -^.^it) /i} ■ (2-28) 

The new Lagrange multipliers, A-^^{t), are related to A^(t) and Xi{t) by 

A,{t) = A,{t) + (3*{t) Sf^(t) = (3*{t) [E,{t) - /i*(t)] + X,{t). (2.29) 

We emphasize once again that the representations for the relevant statistical operator given 
by Eqs. ( p.l2| ), ( p.l7| ), and ( |2.27] ), are equivalent to each other and differ only in the definition 
of the Lagrange multipliers conjugated to the one-particle distribution function. 
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C. Thermodynamic relations 



It is important to note that the relevant statistical operator allows to extend thermo- 
dynamic relations to nonequilibrium systems (see, e.g., W^)- The key quantities are the 
Massieu-Planck function 

m = lnZ,,i(t) = InTr exp |- J] A„(t) P„j (2.30) 

and the nonequilibrium entropy 

S{t) = -Tr {^prei(t) ln^?,eiW} = m + Y: Kit) Pn.it), (2.31) 

m 

which play the role of thermodynamic potentials in the variables {A^(t)} and {P^{t)}, 
respectively. In the case under consideration, using the above given representations for the 
partition function, we obtain formally different but equivalent thermodynamic relations. 
Taking, for instance, the partition function in the form (|2.13|) , we see that the Massieu- 
Planck function can be interpreted as a nonequilibrium thermodynamic potential in the 
variables (3* it), /i*(t), and Xiit). Calculation of its variation gives 

S^t) = -[£- fi*it)N] 6(3*it) + (3*it)N 5/i*(t) - ^ f,it) 5X,it). (2.32) 

1 

In the same representation, the entropy (|2.31|) can be written as 

Sit) = m + [3*{t) [S - /x*(t)iV] + Y: Ai(t)/i(t), (2.33) 



where the self-consistency conditions ( p.l4| ) have been used. Varying both sides of Eq. (p.33|) 
and recaUing Eq. (|2.32| ), we find 



6 Sit) = /3*it) [6£ - /i*(t) 5N] + J2 Ai(t) 6f,it) 

1 

= /3*(t) 68 + [Ai(t) - f3*it) /i*(t)] 6m. (2.34) 

1 

In the last line we have taken into account that 6N = J2i^fi- The entropy may thus be 
considered as a nonequilibrium thermodynamic potential in the variables S and /i(t) or, 
what is the same — in the variables S^^^{t) and /i(t). /,From Eqs. (|2.34| ) it follows that 
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The second equation confirms the fact that thermodynamics determines only the combina- 
tion of the Lagrange multipher Xi{t) and the quasi-chemical potential. 



D. Evolution equations for thermodynamic quantities 

As already discussed, in the approach presented here the total energy is regarded as an in- 
dependent state parameter in addition to the one-particle function . The evolution equation 
for the total energy is trivial: dS/dt = 0. We shall see, however, that the correlation con- 
tribution to the kinetic equation is related to the time dependence of other thermodynamic 
quantities, such as the quasi-temperature and the interaction energy. We will now show 
that the evolution equations for the thermodynamic quantities of interest can be expressed 
in terms of the collision integral ( |2.6|) . 

1. Energy balance 



We start with equation for the interaction energy, Eq. ( 2.15 ). Since the total energy is 
conserved and 

S = {Hy = J2e^f^it)+S,,,{t), (2.36) 



we immediately obtain the balance equation 



dt -T^^ihit)- (2.37) 

In dealing with nonequilibrium many-particle correlations, it is convenient to introduce the 
correlation energy, £^corr(^); which is defined as 

^corr(i) = Ant(i) - ^HfI^) = ^ " E ^1 /l(^) " ^HfI^), (2-38) 

1 

where 
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^hfW = jE(12|V|12)ex/iW/2W = ^Esf'W/iW (2.39) 



is the Hartree-Fock contribution to the total energy. /^From Eqs. ( ^.37| ) and ( |2.38| ) follows 
the balance equation 

'^^corr i^) 



dt 



-E^iWAW (2.40) 

1 

which differs from Eq. ( p. 37] ) in that now the single-particle are given by Eq. (|2.24|) , i.e., 
they involve the exchange Hartree-Fock term. 

2. Equation for the quasi-temperature 

In principle, the evolution equation for (3*{t) may be derived from the equation of state, 
P*(t) = P* {£, {f{t)}), where the second argument indicates that P* is a functional of /i(t). 
This way, however, is not appropriate because we have to calculate the functional derivative 
of f3*{t) with respect to the one-particle distribution function. It is more convenient to make 
use of the self-consistency conditions (|2.19| ) by differentiating them with respect to time. 
This gives 

Tr|^/,UAW, Tr|^J?'U^?£|W. (2.41) 



Since the relevant statistical operator, when taken in the form (|2.17|) , depends on time 
through the Lagrange multipliers /3*(t) and A^(t), we may write 

dt 6p*{t) dt Y ^^i(^) 9t ' ^ ' ' 

With the aid of Eqs. ( p. 17] ) and ( |2.18| ), the variations of g^^ii^) can easily be calculated and 



then Eqs. ( p.41| ) are transformed to (for brevity, the argument t is omitted) 

(A.^')f /.)^^ -A. (2.43, 

{H',H')'-^.pH'jy-ll.-^, (2.44, 
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where we have introduced the correlation function of two dynamical variables: 

(i, b) = £ dx ( Ai g^^, AB g;j^)^^^ (2.45) 
with AA = A — (^)i.er Equation ( |2.43| ) can formally be solved for dA-^^/dt to give 



1 

where is inverse to the correlation matrix 



^^-Exr;|(A..i/')f (2,46) 



Xir = (/i,/i')- (2-47) 



Substituting the expression ( |2.46| ) into Eq. ( |2.44| ) and recalling the balance equation ( p.37D , 
we obtain 

^4i:(^,/0xr.;A. (2.48, 

with the following notation: 

C = {h\ H') - Y: {h', A) Xu' (/i', H') . (2.49) 
11' 

It is easy to check that C in invariant under transformations H' H' + J2i Q^i/i with 
arbitrary coefficients aj^. This property allows one, for instance, to replace H' in Eq. (p.49|) 
by the operator n'{t) [cf. Eqs. ( ^1231) ]. 



3. The entropy production 



We finally derive the entropy balance equation. Recalling Eq. ( p.34| ), we write 



The first term on the right-hand side is zero since the total energy is conserved. Taking also 
into account that df\/dt = and J2i A = dN/dt = 0, we arrive at the equation 

^ = EAiW/i(t) (2.50) 
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which determines the entropy production in the system. We have already noted that, in 
thermal equilibrium, the Lagrange multipliers Xi{t) are equal to zero. In addition, the 
collision integral is also zero in thermal equilibrium. Consequently, the entropy production 
given by Eq. ( |2.50| ) is at least of second order in the deviations from equilibrium, as it should 
be. 



III. THE NONEQUILIBRIUM STATISTICAL OPERATOR 

To proceed beyond the formal thermodynamic relations, we have to calculate the collision 
integral ( p.6|) . In other words, we have to find a solution of the von Neumann equation ( p.4|) 
for the statistical operator g{t) in terms of the state parameters. Since Eq. (|2.4]) is a dif- 
ferential equation with respect to time, one has to give the statistical operator at some 
initial time or require the g{t) to satisfy some boundary condition for t — oo. An 
appropriate choice of the initial or boundary condition depends on the physical situation 
under consideration. Formally, one may consider the initial statistical operator g{to) or the 
limiting statistical operator in a distant past to be arbitrary. We will, however, use the 
initial condition f?(to) = f?rei(^o) some time and the boundary condition that the true 
statistical operator coincides with the relevant statistical operator in a distant past. The 
choice of these special conditions, however, is taken only for simplicity since it is of little 
consequence for the long-time behavior of the kinetic equation, which is of interest to us 
here. 

Assuming that 

^(^o) = QvcM (3-1) 
and using the decomposition (p.22|) of the Hamiltonian, it can easily be verified that the von 
Neumann equation ( |2.4| ) is equivalent to the integral equation 



to 



dt'U,it,t') 4 \git'),n'it')] f/o^(t,t'), (3.2) 



ih 
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where 



(3.3) 



is the unperturbed evolution operator; the symbol exp^ {■ ■ ■} stands for the time-ordered 
exponent. Equation (|3.2|) is still exact. If ?i'(t) is treated as a small perturbation, an ap- 
proximate solution of Eq. ( |3.2| ) can be found by an iterative procedure. We shall restrict our 
discussion to the second-order non-Markovian Born approximation in the collision integral, 
which leads, in the standard approach 0, to the Levinson kinetic equation. 

To calculate the second-order collision integral ( p.6| ), we need the statistical operator 
g{t), correct to first order in the perturbation Ti.'. We shall show later than the time 
derivative dg^^iit') / dt' in Eq. ( |3.2| ) is at least of second order in 7i' and, consequently, 
can be omitted. We next note that the interaction term enters explicitly into the relevant 
statistical operator [see, e.g., ( |2.27D ]. Therefore the leading interaction contribution to the 
commutator [f?rci(i^')) '^o('^')] is at least linear in Ti'] this term must be retained. Finally, 
since the last term in Eq. ( p.2|) is already of first order in the interaction, we may replace 
git') in this term by g^^-^it'). Thus, the nonequilibrium statistical operator, correct to first 
order in 7i', is given by 



Qit) = Qroxit) - / dt'U,{t,t') - U^,{t'),n,{t')] Ul{t,t') 
Jto in L J 

- f dt'u,{t,t')\^W^,{t'),n'{t')\ ul{t,t'). 



(3.4) 



ito " ih 

Another solution of the von Neumann equation can be obtained by imposing the bound- 
ary condition that the true nonequilibrium statistical operator coincides with the relevant 
statistical operator in a distant past. This solution follows easily by adding to Eq. (2^) an 
infinitesimally small source 



dgjt) 
dt 



1 

ih 



g{t),noit)+n'{t)] = -e{g{t) - g,,,{t)} 



(3.5) 



where e — +0 after the calculation of averages with g{t). It is important to note that the 
inclusion of the source term into the von Neumann equation does not violate the energy 
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conservation. This can be seen by taking the product of Eq. (|3.5|) with H and calculating 
the trace. Then we obtain 



d£{t) 
dt 



(3.6) 



The right-hand side of this equation is zero due to the self-consistency condition for the total 
energy in the relevant ensemble. 

Analogously to the derivation of Eq. ( |3.4| ), we use Eq. ( |3.5| ) to derive the first-order 
statistical operator in the form 



-oo III '- -' 

* rft'e-(*-*')L/o(t,t')^[&ci(0,^'(0] U^oit,t'). 

-oo III ^ J 



(3.7) 



This expression can also be interpreted as a rule for passing to the limit ^ —oo in Eq. 
since the factor exp{— £:(t — t')} provides the regularization of the integral. 



IV. THE NON-MARKOVIAN COLLISION INTEGRAL 



We now turn to the calculation of the collision integral (|2.6| ) using the explicit expres- 
sion ( p.4|) for the statistical operator. First we will show that the term Qj.f,\{t) in Eq. (|3.4| ) does 
not contribute to the collision integral. Note that the obvious identity ([/i, In f)j,gi(t)])*gi = 
and Eq. ( p. 17] ) give 



Since = 0, we find that ([/i, -f^'])*ei ~ 0- Thus the collision integral ( |2.6| ) is at 



least of second order in the interaction. This allows us to show that the time derivative 
dg{t')/dt' in Eq. ( |3.2| ) is also of second order in H' and, consequently, it does not contribute 
to the first-order expressions ( |3.4| ) and ( |3.7| ). We use the fact that the relevant statistical 



operator depends on time through the state parameters or through the conjugated Lagrange 
multipliers. For instance, we may assume that Qj-giit') = Qreii^ ^ {f {'t')}) ■ Then, since the 
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total energy is conserved, 

dV ^ 5 hit') dt' V ^fiit') 
whence it follows that dg^^t') / dt' is at least of second order in the perturbation, as was to 
be proved. 

The last two terms in the expression (p.4|), when substituted into Eq. (p.6[), lead to the 
decomposition of the collision integral 

A(t) = /f(t) + jf(t), (4.1) 

where 

I^{t) = -l£dt'Tr{[f/t(t,t')[/i,i^']f/o(i,0,^'(0]f^rei(0}, (4.2) 

I?it) = ^£^/t'Tr{f/t(t,t')[/i,i^1f/o(^,0[^rci(0,^o(0]}- (4.3) 

We shall see later that the term /f (t) is nothing but the Levinson collision integral. The new 
term, /p(t), is due to many-particle correlations in the ensemble described by the relevant 
statistical operator. If the one-particle distribution function fi{t) is taken as the only state 
parameter, the relevant statistical operator ( p.l7| ) does not involve the term with H'; hence, 
Qreii^') commutes with ii-Qit') and /f (t) = 0. 



A. The collision contribution 



We now calculate the term /f (t), Eq. ( |4.2| ), in the non-Markovian Born approxima- 
tion. Since this term is explicitly of second order in the interaction, the relevant statistical 
operator, Eq. ( |2.27| ), can be approximated by 



&(^) = exp |-^Ai(t) ajaij j Tr exp |-^Ai(t) a\a^ . 



(4.4) 



Formally, this statistical operator describes a nonequilibrium ideal quantum gas, so that the 
self-consistency condition for the Lagrange multipliers A^(t) reads 



m = Tr {a\a,gUt)} = r/, , (4.5) 

whence it follows that 

A.(.^.(iif)). 

The time dependence of the operators in Eq. (|4.2|) can be worked out by using the following 
properties of the evolution operator ( pl3|) : 

uUt,t') Uo{t,t') = e-^"i ai, f/o^(t,t') a\ U^{t,t') = e*^i 4, (4.7) 

where 

uj,{t, t') = \f dt" E,{t") = ^-l{t-t') + \ f dt" (t"). (4.8) 
n Jt' a h Jt' 



Taking into account that the statistical operator (|4.4| ) admits Wick's decomposition of the 
averages, a simple algebra gives 

^tit) = -i E I(12|V^|1'2'),J' fdt' cos [Auj,,^Mt,t')] ^i2,i'2' {{fim , (4.9) 

where 

Au,^^,,^,{t,t') =u,{t,t')+u^{t,t')-uj,,{t,t')-uj^,{t,t'), (4.10) 

and the functional JF^2 i'2' ({'Z'}) is defined for any set of single-particle functions (f^ as 

•^12,l'2'({^}) = </'lV'2^1'^2'-<^l<^2</'l'</'2', <^l = l±V'l- (4.11) 

In the context of Eq. (|4.ij|), the functional i'2' ({/(^')}) is nothing but the gain- loss term 
which appears in quantum collision integrals. Expression ( [4.9|) only differs from the original 
Levinson collision integral in that the quantities ^^^(t, t') in Eq. ( [4. 10|) involve the exchange 
Hartree-Fock term. 
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B. The correlation contribution 



To calculate the second-order correlation contribution to the collision integral, Eq. ( [4. 3D , 
we expand the relevant statistical operator ( |2.27|) in Ti', keeping only the first-order terms. 
This gives 

Oroiit) = [l-f3*{t) I'dx [gUt)Y {n'it) - {n%) [gUt)Y^] oUt)- (4.12) 

Here the symbol {'H')q stands for the average with the statistical operator ( [4.4|) . Having the 
above expression, we calculate the commutator appearing in Eq. ([4.3|) (for brevity the time 
argument t' is omitted): 



121'2' 



where 



A^i2,i'2'(i) = E,{t) + E,{t) - Ey{t) - E,,it), 



(4.13) 



(4.14) 



AAi2,i'2'(^) = Ai(t) + A^it) - A,,{t) - A^,{t). (4.15) 

We next substitute Eq. ( [4.13|) into Eq. ([4.3|) and use Wick's theorem to calculate the aver- 
age. In the final result it is convenient to eliminate AA^2i'2'- To this end, we introduce a 
functional 



^12,1'2' ({^}) 

V^l V2 VV ¥^2' 

Then, recalling Eq. (|4.6| ), it can easily be verified that 

AAi2V2'(i) = -ln/Ci2V2'({/(^)}) 



(4.16) 



(4.17) 



Omitting a simple algebra, we present the final expression for the correlation term in the 
collision integral: 
1 

21'2' 



^2 E K12|V|1'2'),J^ / dt' COS Au,,^Mt,t' 



to 



X ^^"^^^TKT^ •^i2,i'2' ({/(^ )}) • 



(4.18) 



in/Ci2,i,2' {{fim 

It is similar to the collision term ( [4.9|) but contains the additional factor in the integrand 
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C. The full collision integral and its properties 



Due to the similarity in structure, the two contributions, Eqs. ( [4.9|) and (^4.18|) , are 
conveniently combined into a single expression 



hit) = -i E Kl2|l^|l'2')ej' fdt' cos [Au,,^Mt,t'] 

ft 01/0/ in 



21'2 



||y}-^i2,i'2'({/(0}), (4.19) 



In/Ci2,i'2' {{fit 

which can be written in a more elegant form by using the relation 

f3*{t) AE,, ,,,,{t) = - InlC,,^,,,, {{F{t)}) , (4.20) 

where 

^'^^^ ^ exp{(3*{t) [Ei(t)-^*(t)]}±l ^^-^^^ 



is the one-particle distribution function in the quasi-equilibrium ensemble [cf. Eq. (|2.21| )]. 
With Eq. ([4.20|) , the expression ( [4.19|) becomes 



hit) = E I(12|^|1'2'),J' fdt' cos [Au;,,^Mt,t' 



This collision integral has some remarkable properties. First, it vanishes in the quasi- 
equilibrium state where fi{t) = Fi{t). In this case the collision contribution and the correla- 
tion contribution cancel each other. In particular, = in complete equilibrium since {t) 
goes over to the equilibrium distribution function as l3*{t) ^ 1/T and fi*(t) fi, where T 
and fi are the equilibrium temperature and the equilibrium chemical potential, respectively. 
It should be emphasized that the collision term ( [4.9| ) alone does not vanish in thermal equi- 
librium, which is the grave disadvantage of the Levinson-type kinetic equations. Another 
important property of the collision integral is its asymptotic behavior as t — 00. This 
stage of the evolution can be described on a large time scale, so that we may pass to the 
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Markovian limit. To analyze this case, it is convenient to return to the expression ( [4.191) . 
Putting there f{t') ^ f{t), F{t') F{t), and then performing the limit — oo with the 

regularization factor exp {—6{t — t')}, we find that the correlation contribution vanishes due 
to the fact that now the cosine term is replaced by the delta function 6 (^AEi2 i'2'(^)/^)- 
a result, we get the well-known Uehling-Uhlenbeck collision integral 

m = -^Y:J{12\V\1'2')J 6 (^^^1^^ (^fjj^j^, _ A/,/,/,)^, (4.23) 

where we have used the definition of ^i2i'2') Eq. (|4.11| ). We would like to draw attention 
to the role of the correlation term in the collision integral. Although this term goes to 
zero in the long-time limit, the non-Markovian expression ( [4.19|) is constructed such that 
the interplay between collisions and correlations is precisely the reason why the Markovian 
regime arises. It should be noted, however, that beyond Born approximation, for instance, 
in the T-approximation for the collision integral, the correlation term does not go to zero in 



the Markovian limit [O]. 



D. A simplified version of the non-Markovian collision integral 



Because of the presence of the /C-functional, the full collision integral (|4.19|) has a more 
complicated structure than the Levinson term ( [4.9|) . Having in mind practical applications 
of the scheme developed here, it makes sense to formulate a simplified version of the non- 
Markovian collision integral which, nevertheless, retains the main properties of the full 
expression ( [4.19| ). Let us approximate the correlation term ( [4.18| ) by its value in the quasi- 
equilibrium state described by the statistical operator ( |2.21| ) . In the case of weak interaction, 
this approximation means that we put f{t') F{t'). Then, recalling Eq. (|4.20|) , we obtain 



n 2i'2' 



to 



dt' COS Au;i2,l'2'(^,^') •^12,1'2' ({^(0}) • 



(4.24) 



Now the expression ( |4.19 ) takes a simpler form 
1 



21'2' 



J2 |(12|y|l'2'),j' / dt' cos Auj,,^Mt,t' 



to 
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X 



-^i2,i'2' ({fit')}) - -^i2,i'2' i{F{t')})] ■ (4.25) 



Obviously this collision integral vanishes in the quasi-equilibrium state and, consequently, in 
complete equilibrium. Another important point is that, in the Markovian limit, Eq. ( [4.25[ ) 



reduces to the Uehling-Uhlenbeck collision integral ( [4. 231) since, in this limit, the contri 



bution from JF^2 i'2' ({-^(^')}) vanishes, as can easily be verified with the aid of Eqs. ([4.11|) 



and (|4.21 ). Recently, an "improved" version of the Levinson collision integral was proposed 



in the Green's function method on the basis of approximate solution of a Dyson equation 
with initial correlations 0. It differs from Eq. ( 4.25 ) in that the single-particle energies 



did not involve the Hartree-Fock corrections and the correlation term was approximated by 
jF^2,i'2' ({Z'"*^''''})) where fi"^^ is the distribution function in complete equilibrium. Physically, 
the drug of choice for a simplified non-Markovian collision integral is the expression (f4.25 ) 



which involves effects of running correlations, whereas replacing the quasi-equilibrium dis- 
tribution function F^{t') by fi"^^ implies that the state of the system is close to complete 
equilibrium. 

V. BALANCE EQUATIONS IN THE NON-MARKOVIAN REGIME 

Having the explicit expression ([4.19|) for the collision integral, it is of interest to analyze 
in more detail the balance equations derived in Section II. 

A. Energy balance. The nonequilibrium correlation energy 

We have already shown that the non-Markovian kinetic equation with the collision in- 
tegral ( |4.19| ) has the correct equilibrium solution. Now we want to demonstrate that the 
right-hand side of Eq. ( |2.4CI| ) can be represented as a time derivative, i.e., the non-Markovian 
kinetic equation is consistent with the energy conservation. 

First we will prove that any collision integral of the form 

rt 
I to 



hit) = -Y. f dt' cos [A^i2,1'2'(^,^')] Gi2,i,2'(i') (5.1) 
21'2' *0 
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conserves the total energy, if the function i'2' (^) has the symmetry properties 

^12,l'2'(^) ~ ^21,1'2'(^) ~ ^12,2'1'(^)' ^12,1'2'('^) ~ ~^1'2',12(^) (5-2) 



and the function Alj^^j i'2'(^5^') satisfies the conditions 
d 

^ ^^^12,1'2' t') = ^^12,1'2' (^) ) ^^^12,1'2' ^) = 0, 



(5.3) 



where Ai?^2i'2'(^) given by Eq. (|4.14| ). The proof is as follows. Multiplying Eq. (|5.1| ) by 
E^^it) and then summing over the quantum numbers 1, we obtain 

Y,E,{t)I,{t) = -- ^E.^Mt) fdt' cos[Aa;i2,i,2,(t,0]G'i2,i,2'(0 

1 121'2' *0 

= ^ (-^ E Hdt' sin [A^i2,i,2,(t,tO] Gi2,i,2'(i')) , (5-4) 

where use has been made of Eqs. ( p.2| ) and ( |5.3| ). Comparison of Eqs. ( |2.4CI| ) and ( p.4| ) shows 
that the non-Markovian collision integral indeed conserves the total energy. As a by- 
product of the proof, we have the following expression for the nonequilibrium correlation 
energy: 

^ /■* r -1 

^corrW = ^corr(i^o) + T E / dt' sin[/\uJ^^ y^,{t,t')\G^2 y^,{t') , (5.5) 

121'2' " 

where i^corr(^o) initial value of the correlation energy at t = t^. 

Turning now to the collision integral ( [4. 191) and recalling the definition of the functionals 
JF and /C, it is easy to check that the symmetry conditions (|5.2|) are satisfied. Thus, the 
collision integral ([4.19|) conserves the total energy. In our case Eq. (|5.5|) reads 



^corr(^) ~ ^corr(^o) + ^^corA'^) + ^^corrl*^)) 



(5.6) 



where 



4 y Ki2|v|i'2'), 



AC„W = Tr E l(12|l^|l'2') 

"ill -^21/2' 



1 2 
ox I 



/ dt' sin 


^'^12,1'2'(^' 


'to 






1 dt' sin 


^'^12,1'2'(^) ^0 


'to 



/3*(tOAE,2,r2^(tO ^ 
\r.lC,,,,A{f{t')}y '''''' ^^^^^^^ 



(5.7) 



(5.8) 
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We have separated the time-dependent contribution to the correlation energy into two parts 
which have different physical interpretation. The term ( |5.7| ) can be regarded as the collision 
contribution to the correlation energy. An analogous term was derived previously from the 
Levinson kinetic equation [|15|,|T6| and by the Green's function method ||^. The term ( p.8| ) 
arises due to collective (correlation) effects. In thermal equilibrium these two terms cancel 
each other so that the correlation energy does not depend on time. 

It is interesting to note that the simplified non-Markovian collision integral (|4.25| ) has 
the form (|5J] ), where Gi2,i'2' is the difference of two jF-functionals, each of which satisfies 
the symmetry conditions ( |5.2|) . We may thus conclude that the collision integral (|4.25|) 
conserves the total energy. The approximate correlation energy is given by 
1 ^ ,0 /•* 



121'2' 



^corrW = ^ \{12\V\1'2')J J^^^ dt' Sin [^UJ,,^y,,{t,t') 

X {-^i2,i'2' iifim - -^i2,i'2' am})} ■ (5.9) 



Exactly the same expression follows from Eq. ( |5.6| ) if the term ( p.8[ ) is replaced by its value 
in the quasi-equilibrium state. 



B. Non-Markovian equation for the quasi-temperature 



In general, the quasi-temperature evolves in time according to Eq. ( |2.48|) . Within the 
framework of non-Markovian Born approximation, the correlation function {H, f\) in this 
equation can be replaced by {HQ{t), f^) since the collision integral is already of second order 
in the interaction. Then a little algebra shows that Eq. ( |2.48| ) reduces to 

dl3*{t) 1 



With Eq. (^4.19 ), this is written in the expanded form as 
d^*{t) 1 



dt 



Afi^ C{t) 



121'2' 



X / dt' COS AuJ^2 1'2'i^^'^') 
'to L 



P*{t')AE,,^Mt') 
ln/C,2V2'({/(i')}) 



(5.10) 



•^i2,i'2'({/(0})- (5.11) 
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The correlation function C{t) is given by Eq. ( p.49|) and can be calculated, in the leading 
approximation, by using the statistical operator ( [4 .41 ) which admits Wick's decomposition of 
averages. The Lagrange multipliers A^^ (t) can then be expressed in terms of the one-particle 
distribution function by means of Eq. (f4.6|). After some algebra, one obtains 



c{t) = \y: mivmef mmfvit)fAt) ^'''P' ^^^amu^ ^^-^2) 

^ 121'2' ll/l'-Jj J 

whence it follows that C{t) > 0. Now the kinetic equation ( |2.5|) , together with the ex- 
pression ([4.19 ) for the collision integral and the evolution equation ( 5.1 1| ) for the quasi- 
temperature, form a closed set of equations describing non-Markovian relaxation processes 
in the system. 



C. Entropy production in the non-Markovian regime 



Of special physical interest is the entropy equation (p.50|) . To second order in the interac- 
tion, the Lagrange multipliers Xi{t) can be expressed in terms of the one-particle distribution 
functions, f^{t) and F^(t), with the aid of Eqs. ( p?29D , (TO, and (|]21]). Eliminating A^ and 

P*{E-j^ — fi*), we obtain 



Ai(t) = In 



(5.13) 



Substituting this expression, together with Eq. ([4.22|) , into Eq. (|2.5CI|) and then making use 
of the symmetry of the collision integral under permutations of the single-particle quantum 
numbers, the entropy production takes the form 



dS{t) 
dt 



Y: Kl2|V^|l'2'),jMn 

t |- 

X / dt' COS A(x';^2 l'2'(^) 
'to 



N ln/Ci2,i.2-({F(tO}) ] ^ 

i^-in^.,..({/(to})r--^^^^^^»- ^'-''^ 



Although the right-hand side of his equation involves memory effects, the entropy production 
is identically zero in thermal equilibrium, as it must be. To understand the behavior of the 
entropy production in the Markovian limit, we put f{t') ^ f{t), F{t') ^ F{t), and then pass 
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to the limit ~^ "Oo inserting the factor exp{—e{t — t')}. As we have already noted, such 



a procedure leads to the appearance of the delta function 5(Ai?^2 i'2'(^)/^) place of the 
oscillating cosine term. Due to Eq. (|4.2CI| ), we may then put /C12 i'2' ({-^(^)}) = 1- Finally, 
it is convenient to eliminate the jF-functional be means of the relation 



^12,1'2' ({/}) = /l/2/l'/2' (/Ci2,l'2' ({/}) - 1) 



(5.15) 



which follows from Eqs. ( [4.11| ) and (f4.16|) . As a result of these manipulations, we find the 
entropy production in the Markovian limit: 



121'2' 



h 



X 



ium) - 1 in/Ci,i,2' ium) /i/2/i'/2' 



(5.16) 



It depends only on the one-particle distribution function and is positive, since (a; — 1) Inx > 
for X > 0. The result ( |5.16|) agrees with the well-known expression for the entropy production 
in a weakly interacting quantum system described by a Markovian kinetic equation (see, 

e.g., m- 



VI. GENERALIZATION TO SPATIALLY NON-HOMOGENEOUS SYSTEMS 

Here we shall briefly touch upon the extension of the foregoing treatment to spatially 
non-homogeneous systems. In such cases, many-particle correlations associated with the 
energy conservation can be incorporated by taking mean energy density S{r,t) = {H{r)Y 
as a state variable, together with the one-particle Wigner function /^(r, p, t) = (/^.(r, p))*. 
The energy density operator, H{r), is defined through the relation 

H = J drH{r), (6.1) 

and the operator corresponding to the Wigner function is given by 

/Jr,p) = J cixe-P-^/'^V^t(r_x/2)^Jr + x/2), (6.2) 
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where ipai'^) and ipa^r) are the second-quantized field operators; the Greek indices a = {i, a) 
specify the type of particles {i) and the spin state (a). Now the relevant statistical operator 
takes the form [cf. Eq. (g]T2|)] 
1 



^rel(^) 



exp <^ - / rfr/3*(r,t) 



-E 



i 

dp 



A„(r,p,t)/^(r,p) L 



(6.3) 



(27r^)3 ^ 

where n^ir) = X^o- V^i(r)'?Aa(r) are the particle-number density operators for the species. The 
local inverse quasi-temperature /3*(r,t) and the Lagrange multipliers X^{r,p,t) are to be 
determined from the local equations of state 



£{r, t) = Tr {^(r) g,,,{t)} , /Jr, p) = Tr {/Jr, p) g,,,{t)} . 



(6.4) 



Just as in the spatially homogeneous case, the local quasi-chemical potentials /i*(r,t) in 
Eq. ( |6.3| ) can be eliminated by re-defining the Lagrange multipliers X^{r,p,t). Another 
possible way is to define the /i*(r, t) in the local-equilibrium state described by the statistical 
operator 

1 



exp<^ - / dv[3*{v,t) 



i7(r)-5:/i*(r,t)n,(r) 



(6.5) 



Then the quasi-chemical potentials are to be determined from the local equations of state 



ni(r,t) = Tr{ni(r) ftocW) 



(6.6) 



where n-{r,t) = (n,-(r))* are the mean particle-number densities for the species. The latter 
definition of the quasi-chemical potentials is typical for the hydrodynamic description of 
transport processes. 

In the spatially non-homogeneous case, the basic evolution equations are the kinetic 
equation for the Wigner function 



|/Jr,p,t) = i([/Jr,p),^])* 
and the local conservation law for the mean energy density 



(6.7) 
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|f (r,t) = 1 {[H{r),H]y = -V ■ (J(r))*, (6.8) 



where J(r) is the energy flux operator. To calculate the averages appearing in Eqs. 
and ( |6.8|) , the nonequilibrium statistical operator has to be found as a functional of the 
state variables. This can be done by solving Eq. ( p.2|) in some approximation with the 
relevant statistical operator given by Eq. ( |6.3| ). If, for instance, the interaction term in the 
Hamiltonian can be regarded as a small perturbation, then Eq. ( p.2| ) can be solved by an 
iterative method which is similar to the procedure used in Section III. There are, however, 
some new features which are specific to spatially non- homogeneous systems. First, now the 
average on the right-hand side of Eq. ( |6.7| ), when calculated with the relevant statistical 
operator, is not zero and the kinetic equation for the Wigner function has the form 

^ /Jr, p, t) + D Jr, p, t) = 4(r, p, t), (6.9) 

where 

D^ir,p,t) = ^{[ar,p),H])l, (6.10) 



is a generalized drift term which involves the correlation contributions [0]. The collision 
integrals /Q,(r, p, t) are determined by the integral terms in Eq. (|3.2|) . The second important 
point is that, in the spatially non-homogeneous case, the right-hand side of Eq. ( p.2| ) contains 
terms with gradients of the Lagrange multipliers, /5*(r, t), /i*(r, t), and A*(r, p, t). Finally, 
the collision integrals /^^(r, p, t) are, in general, non-local functionals of the Wigner function. 
If the state parameters vary little over the range of the interaction potential and the mean 
de Broglie wavelength, which is typical for real situations, then the non-locality effects can 
be incorporated by using an expansion of averages in Eqs. (|6l7| ) and (|6l8| ) in powers of 
spatial gradients. Within this scheme, the evolution equations derived in this paper may be 
interpreted as transport equations in the local approximation, where all the gradient terms 
are neglected. 
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VII. CONCLUSIONS AND OUTLOOK 



We now summarize the main implications of the above analysis. The first point is that a 
non-Markovian kinetic equation conserves the total energy and has the correct equilibrium 
solution only if the collision and correlation effects are incorporated in a self-consistent 
way, for instance, within the same approximation in the interaction between particles. We 
have also seen that it is precisely the interplay between correlations and collisions that is 
responsible for the long-time behavior of the non-Markovian collision integral. 

Another important feature of the outlined approach is that the kinetic equation for the 
one-particle distribution function is supplemented by the equation for the quasi-temperature, 
Eq. (15.101 ), which describes the "slow" evolution of the system. This representation for non- 



Markovian kinetics differs from standard schemes, say the Green's function method |18-20 



where the goal is to derive a closed kinetic equation, so that the collision and correlation 
effects have to be described in terms of single-particle quantities. The inclusion of long- 
lived many-particle correlations into the Green's function formalism is a rather difficult 



problem and there are only first steps in this way pT|-p3|. On the other hand, the Green's 
function technique provides a powerful tool for calculating the collision contribution to a 
kinetic equation beyond the Born approximation. Thus an interesting point would be a 
unification of the Green's function method and the density operator method to develop the 
self-consistent non-Markovian quantum kinetic theory involving the correlation effects and 
the quasiparticle damping. 

It is significant that the approach outlined in this paper is non-perturbative in external 
fields. If the external field does not directly affect interactions between particles, then 
the nonequilibrium statistical operator will have the form ( p.4|) , but the evolution operator 
f/g(t, t') will now involve the field effects which can be taken into account exactly (not in 
terms of perturbation theory). It should also be noted that, in the presence of an external 
field, the energy of the system is not conserved and the trivial equation dS/dt = has 
therefore to be replaced by a balance equation including the work produced by the field. 
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